Continuously decoupling single-qubit operations from a perturbing thermal bath of 

scalar bosons 
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We investigate the use of continuously-applied external fields to maximize the fidelity of quan- 
tum logic operations performed on a decohering qubit. Assuming a known error operator and an 
environment represented by a scalar boson field at a finite temperature, we show how decoherence 
during logical operations can be efficiently reduced by applying a superposition of two external 
vector fields; one rotating orthogonally to the direction of the other, which remains static. The 
required field directions, frequency of rotation and amplitudes to decouple noise dynamically are 
determined by the coupling constants and the desired logical operation. We illustrate these findings 
numerically for a Ifadamard quantum gate and an environment with ohmic spectral density. 
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A quantum computer, when finally built, will be more 
efficient than current classical computers to solve certain 
kinds of problems The theory of quantum informa- 
tion processing generally takes advantage of the inherent 
parallehsm exhibited by unitary operations on quantum- 
state superpositions. The terms of these linear combina- 
tions are tensor products of quantum bits, or "qubits" 
0, which, linearly superposed, result in states with the 
desired properties of entanglement and interference 
In principle, the choice of an appropriate external field 
would guarantee a correct dynamics for the system, se- 
lected among those exhibiting unitary symmetry. How- 
ever, during the actual quantum evolution of the sys- 
tem, since it cannot be completely separated from its 
environment, the unitary symmetry breaks down. The 
consequent decay of the quantum state purity is a man- 
ifestation of the ubiquitous phenomenon of decoherence 

i- 

There are at least three major classes of strategic de- 
vices proposed to counteract the deleterious and unavoid- 
able effects of decoherence: quantum error correcting 
codes , decoherence- free subspaces and subsystems Q , 
and dynamical decoupling 0, H B ■ Because the first 
two of these strategies require more than one physical 
qubit to protect each logical qubit, dynamical decoupling 
is the simplest of the three, since it requires, in princi- 
ple, only controllable external fields to directly protect 
each physical qubit. Even without precise knowledge of 
the error structure and strengths, the pulsed dynamical- 
decoupling scheme is effective, but usually employ an 
articulate time sequence of external-field pulses which, 
for experimental implementations, requires sophisticated 
control procedures. Moreover, the pulses have to be so 
short as to start and finish well within the environmental 
correlation time interval, so the field intensities involved 
must be high. Initial attempts to use continuously- 
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app lied fields instead of pulses have appeared recently 
and these preliminary analyses show that, although 
pulses are not necessary for protecting against the effects 
of particular error structures assumed known, employ- 
ing fast control cycles is inevitable. From the practical 
point of view of experimental realization, however, the 
successful introduction of continuous control fields into 
the general dynamical-decoupling framework will be rel- 
evant even in the high-frequency regime. 

In this paper we show that, if the error structure and 
strengths are known, by continuously applying a suit- 
able superposition of external fields it is possible to re- 
alize, with high fidelity, any logical operation on a qubit 
weakly perturbed by an environment represented by a 
scalar boson field at a finite temperature. The general 
formulation of the dynamical decoupling method [1, [l^l 
gives a clear geometrical interpretation of the error pro- 
tection when explicitly applied to the simple case of a sin- 
gle qubit driven by continuous fields. Hence, under these 
circumstances, we have found that decoherence during 
logical operations can be efficiently reduced by applying 
a superposition of two external vector fields: one that 
rotates orthogonally to the direction of the other, which 
is a static vector field. The amplitudes, frequency, and 
directions of these fields are determined by the intended 
quantum logic operation, the error structure, and the 
characteristics of the environment. 

For the purpose of starting an investigation on the use 
of continuous control fields to protect against errors dur- 
ing qubit operations, we make the simplifying assump- 
tion that the interaction between the qubit and its en- 
vironment is sufficiently weak that perturbation theory 
is applicable. For arbitrary coupling strengths to un- 
desirable but known terms in the qubit-system Hamil- 
tonians, without taking into account environmental de- 
grees of freedom that are to be traced out, there are the 
optimal-control approaches of Ref. In the present 

work we imagine a situation in which the qubit is far from 
other qubits, but is acted upon by controllable external 
fields and is subject to residual influences from the sur- 
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round environment. For this particular setting, Ref. fl3] 
presents an optimal control theory where the environ- 
mental perturbations are modeled by a Lindblad master 
equation, which is valid in the Markovian limit [isj . Since 
here we treat the details of the controlling cycles occur- 
ring during the correlation time of the fields representing 
the environment, we must, although using perturbation 
theory, derive a non-Markovian master equation. 

We start by noticing that the agents coupling the qubit 
to the environment can be represented by the Pauli ma- 
trices. For the sake of simplicity, we represent the action 
of the environment by a single scalar boson field and as- 
sume an interaction Hamiltonian given by 



ifint = {X-(t)B+{X* -(7)3^, 



(1) 



where A is the error vector whose components are com- 
plex and we take B = X^/tSfcOfc, where gk is a complex 
coupling constant for normal mode k with dimension of 
frequency and ak is the operator that annihilates a bath 
quantum in mode k. Here and in the following we use 
units of h ^ I. 

Instead of starting from the Hamiltonian and deriv- 
ing the propagator, which requires dealing with techni- 
cal issues related to time ordering of interaction oper- 
ators, we begin with the most general unitary evolution 
and obtain the Hamiltonian by differentiation. Any time- 
dependent unitary transformation of a qubit can always 
be expressed as 



U{t) = / cos [a{t)] 



u{t) sin [a{t)] , 



(2) 



where / is the 2x2 identity matrix, a (t) is a function of 
time t, u (t) is a time-dependent unit vector, cr = xct^: + 
yay+zaz, and a^, ctj,, and cr^ are the Pauli matrices. The 
corresponding driving Hamiltonian Hij{t) is obtained by 
differentiating Eq. 



Hu{t)=^'^U^{t) = m 



(3) 



with 



n[t) 



da(t) 
dt 



siT?[a{t)]n(t)x 



i[a(t)] cos[a(i)] 

d\i{t) 

dt ' 



du(t) 
dt 



(4) 



where U^{t) is the Hermitian conjugate of U{t). It is im- 
portant to emphasize that because of the results in Eqs. 
(I2])-(|3]) we avoid any time-ordering manipulations and do 
not need the so-called polaron transformation [l^. 

The time-local second-order master equation describ- 
ing the evolution of the reduced density matrix of the 
qubit, in the interaction picture, is written as (l7| 



dpijt) 
dt 



dt'TiB {[Hiit), [Hi{t'), pBPiit)]]} , (5) 



where Hi(t) is the interaction Hamiltonian in the interac- 
tion picture, namely, Hi{t) = W{t)ul{t)HintUB{t)U{t), 



with Usit) = exp {-iHst), He = J^k^kalat, LOk is the 
frequency of normal mode k of the thermal bath, and 
U{t) is as in Eq. Equation ([5]) is valid in the regime 
in which the strength of the coupling, expressed in fre- 
quency units, multiplied by the correlation time of the 
bath operators is much lesser than unity. 

Above, ps is the initial density matrix of the thermal 
bath: 



Pb = -^exp(-/3i/B), 



(6) 



where Z is the partition function given hy Z = 
TiB [exp(— /Ji/s)]. Here, (3 = l/ZcsT, ks is Boltzmann 
constant, and T is the absolute temperature of the envi- 
ronment. 

From the form of the interaction between the qubit 
and its environment, Eq. ([1]), we obtain 



Hi{t) = Bi{t)A{t) ■ a + B]{t)A*{t) ■ a. 



(7) 



where Bi{t) = ul{t)BUB{t) = J2k9kakexp{-iu!kt) and 
we have defined the time-dependent vector A(i) as 

A{t) = Acos[2Q;(t)] -I- [Axu(t)]sin[2a(i)] 

+ u(t) [u{t) - A] {1 - cos [2a{t)]} . (8) 

Substituting Eq. ([7|) into Eq. we obtain 

3 3 



dp lit) 
dt 



a = l 13=1 
3 3 

Q=l /3=1 



(9) 



where we have defined 



D^pit)^ f dt'TTBMt)pBb0it')] 
Jo 



(10) 



with the bath vector operator h{t) given by h{t) 
Biit)Ait)+B\{t)A*{t). 

Substituting ^ into PU]) gives 



D^pit) = 2Re A*^{t) / dt'Ap{t')Ii{t - t') 



+A*Jt) / dt'Apit')l2it-t'), (11) 
Jo 



where we have defined 



9k\ 



exp (iujkt) 
exp(/?Wi;) - 1' 



^2{t) = \gk\^ exp (iuJkt) . 



(12) 
(13) 



In the limit in which the number of bath normal modes 
per unit frequency becomes infinite, we define its spectral 
density as 



'^('^) = E \9kf^{^ - ^k), 



(14) 
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with uj €E [0; +00) and interpret the summa- 
tions in Eqs. and (|13p as integrals over uj: 
Ti{t) = /g°° (iti;J(a;) exp(icji)/[exp(/3a;) — 1] and l2(i) = 
duiJ{uj) exp{iu}t). 
If we now write pi{t) — 1/2 + r{t) ■ tr, the Bloch vec- 
tor r{t) is real and, from Eqs. ([9]) and (fTTjl . satisfy the 
differential equation 



dr(t) 



= 4Rc{A*{t)x[{2F{t) + G{t))xr{t)]} 
-2Im[A*(i) X G(t)], (15) 



where we have defined F{t) = dt' A{t')Ji{t - t') and 

G{t)= ^ldt'A{t')X2{t-t'). 

The vectors F(i) and G(t) can be interpreted as time 
averages of the vector A(i') weighted by the functions 
Xi{t — t') and l2{t — t'), respectively. Depending on the 
width of the reservoir spectral density, J(w) in Eq. (fT4|) . 
the weight functions, Ti{t — t') and T2{t — t'), determine 
how much of the recent past history of A{t') effectively 
contributes to its time averages, F{t) and G{t). Hence, 
Eq. (jisp is not restricted by a Markovian approximation. 

We aim at finding U{t) such that Eq. which 
refers to the interaction picture, gives, after a certain 
gate-operation time r, r(T) « i"(0). In the Schrodinger 
picture, this result means that, effectively, the dynamics 
described by Eq. (fT5|) is equivalent to the action of U (r), 
as if the qubit were not perturbed by its environment. 
The total Hamiltonian, in the Shrodinger picture, is given 
by H{t) = Hu{t) + Hb + Hint, where Hu{t) is to be 
determined as prescribed by Eq. ^ after we find U{t), 
Hb is the bath Hamiltonian as described just below Eq. 
([5]), and Hint is the interaction Hamiltonian as given by 
Eq. (HD. 

Because we intend to realize a quantum logic operation 
simultaneously to the protection from errors, we split Hjj 
into two terms, Hij{t) = Ho{t) + Hc{t), where Ho{t) is to 
produce the quantum-gate action, while Hc{t) counter- 
acts the perturbing action of the environment. According 
to the prescription well exposed in Ref . [l^] , the unitary 
operator Uc{t) corresponding to the Hamiltonian Hc{t) 

is to be periodic and satisfy J*" dtU^{t)HintUc{t) — 0, 
where tc < t is the period of Uc{t), that is, Uc{t + tc) = 
Uc{t). Here, for convenience, we choose r as an in- 
teger multiple n of tc, that is, r — ntc, such that 
Uc{t) = I. If we choose a constant He = {2TT/tc)uc ■ c, 
then Uc{t) — I cos{2n'Kt / t) — iUc • (Tsin(2n7rt/r). By the 
same procedure to obtain Eqs. ([7]) and ([5]), we find that 
the integral of U^{t)HintUc{t) over one period tc gives 
zero only if we choose Uc orthogonal to the error vec- 
tor A, which is always possible, even if A has complex 
components. 

In the picture obtained using the unitary transfor- 
mation Uc{t) as given above, we choose Ho{t) such 
that UUt)Ho{t)Uc{t) = I cos 00 — i(T ■ Uo sin 6*0, where 
Uq and ^0 are constants to be determined according 
to the intended quantum logic operation as explained 
in the following. Since U^{t)HQ{t)Uc{t) is time inde- 
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FIG. 1; Numerical solutions for the strategy example de- 
scribed in the text to overcome decoherence and dissipation 
during a Hadamard operation. The dot-dashed line represents 
the fidelity when He = 0, while the solid line is the result for 
n = 5 in Eq. (|18|l and Uc, chosen orthogonally to the error 
vector A, as in the text. The dotted and dashed lines repre- 
sent the n = 5 results when Uc is tilted 10° and 30° toward 
the X axis, respectively. In the inset we show the fidelities at 
t = r as functions of temperature and n. 



pendent, the associated unitary evolution operator is 
given by UQ{t) = /cos(6'ot/T) — zuq • cr sin{9ot/T), giv- 
ing Uo{t) — I cos do — icr- Uq sin 9q at the end of the gate 
operation, t = t. Therefore, if we take U{t) as the com- 
posed unitary operator U{t) = Uc{t)Uo{t), then, at t = r, 
U{t) = / cos 00 ^ ic ■ liosin^o, so that the intended op- 
eration determines our choice of Uo and ^o- Comparing 
these conclusions with Eq. ([2]), we obtain: 

cos[Q!(t)] = — Uc • Uq sin(2ri7ri/T) sin(0o/''') 

+ cos(2ri7ri/T)cos(6lo/r), (16) 
u(i) sin[a(t)] — (Uc x Uq) sin(2n7rt/T) sin(^?o/''') 
+Uc sin{2mTt/T) cos{9o/t) 
+Uo cos{2mrt / t) sin{9o / t) . (17) 

The explicit form of the Hamiltonian Hjj is the one 
already given by Eq. Hjj — ^{t) ■ cr, where the 

applied external field is calculated from Eqs. (|16p and 
([TT]) according to the prescription of Eq. 

n{t) = [{2mr/T) + i9o/T)uc-Uo]uc 

+ {9o/t)[uc X (uo X Uc)] cos{2mrt/T) 

+ {9q/t){uc X Uo) sin(2n7ri/r). (18) 



The first term of Eq. is a static field along the 

direction that is perpendicular to the error vector, as 
discussed above, and the other two terms give a rotating 
field perpendicular to the direction of the static field. 

We have performed numerical calculations for the 
strategy described above by solving Eq. I|15p. To il- 
lustrate our typical results for a concrete example, in the 
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present article we assume an ohmic spectral density with 
a cutoff frequency Wc, namely, J{oj) = r]uje'}cp{—uj/oJc), 
where 77 is a dimensionless constant. Hence, the integral 
versions of Eqs. and can be explicitly calcu- 

lated to give Ii{t) = (?7//32)<If(i)(l + l/(/3wc)-iV^), 
and l2{t) = rjuj^/ {\ ~ iLOct)^ ^ where "^^^^ is the first 
polygamma function. Since, as in the case of pulsed dy- 
namical decoupling ^, ^ Q\ , we have found that the ro- 
tating field, to be effective, must rotate at a frequency 
sufficiently higher than lDc, we take lOcT = 27r in the nu- 
merical calculations. One of the worst situations occurs 
when the error vector has complex components, describ- 
ing, besides pure decoherence, also dissipation. Thus, we 
take A = (4x + iy + 2\/2z)/5, since it is dimensionless 
according to Eq. Thus, Uc = (x— V2z)/V3 is a suit- 
able choice. For a Hadamard operation, we can choose 
uo (x + z)/V2 3. We take -q = 1/16, T = 0.25K, and 
T — 10^ ^"s. With these numbers and the initial condi- 
tion (0/(0) = 1/2 -)- cra;/2, when He = the fidelity gives 
J^(t) = Tt[pi{t)pi{0)] « 0.7199. Now, ff we take n = 5 
in Eq. (flSjl . the fidehty becomes J-{t) w 0.9965. Figure 
1 shows these fidelities as functions of time, together with 
two cases in which Uc is not chosen as orthogonal to the 
error vector, but tilted 10° and 30° toward the x axis. 
In these examples we see that the resulting fidelity func- 
tion is only smoothly sensitive to small-angle variations. 
Incidentally, since the condition for the present scheme 
to work depends only on the angle between and A, 
the error protection is also smoothly dependent on the 
error- vector direction. In the inset, we show the final fi- 
delities, calculated ai t — t, as functions of temperature 
and n with Uc as given above, orthogonal to A; we see 
that lower temperatures are better for increasing fidelity. 



but higher-temperature effects can be compensated by 
higher control-field frequencies. 

The initial conditions are specified by r(0), whose mod- 
ulus, for pure initial states, is 1/2. Thus, all the initial 
conditions for pure initial states can be parametrized by 
two angles: (p G [0, 27r) and 9 S [0, tt]. We have par- 
titioned these angle ranges into 200 and 100 regularly- 
spaced intervals, respectively, and have solved Eq. (fT5)) 
for each of these 20, 000 initial conditions. The worst and 
best cases resulted in fidelities of 0.99622 and 0.99987, re- 
spectively, with all the other variables set as for the solid 
line in Fig. 1. 

In summary, we have shown that it is possible, using 
continuously-applied external fields, to realize one-qubit 
high-fidelity quantum logic operations and, simultane- 
ously, protect these operations against environmental er- 
rors. Our proposed continuous dynamical decoupling is 
less vulnerable to errors produced by non-ideal driving 
than the pulsed version: our method requires interven- 
tion only at the initial and final steps of the quantum 
logic operation, while the pulsed method demands so- 
phisticated timing control for each of the various neces- 
sary pulses. 
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